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Abstract. The paper examines one-dimensional total variation flow equation with Dirichlet 
boundary conditions. Thanks to a new concept of "almost classical" solutions we are able to 
determine evolution of facets - flat regions of solutions. A key element of our approach is the 
natural regularity determined by nonlinear elliptic operator, for which x 2 is an irregular function. 
Such a point of view allows us to construct solutions. We apply this idea to implement our 
approach to numerical simulations for typical initial data. Due to the nature of Dirichlet data any 
monotone function is an equilibrium. We prove that each solution reaches such steady state in a 
finite time. 
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1 Introduction 

The equation which is the topic of this paper 



Ut - ^ (sgn(u a )) = 0, u(a) = a b , u(b) = a e . (1.1) 

is a one-dimensional example of the total- variation flow. The motivation to study this problem is 
twofold: a) image analysis, see BROF1 . HACU . Kill : b) crystal growth problems, see HAGH . iTTal . 
[GGK], [Ma|. There are different physically relevant models, where a similar to ours surface 
energy appears, but the corresponding evolutionary problem is not necessarily set up, see e.g. 

Equation (11.11) may be interpreted as a steepest descent of the total variation, i.e. we can write 
(11.11) as a gradient flow u t E —dE{u) for a functional E. This is why we can apply the abstract 
nonlinear semigroup theory of Komura, see llBrl . II Ball , to obtain existence of solutions. This has 
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been performed by ED, HKGl - llGGKl and also by llABCll . llABall . llACDl . llBCNl . However, 
the generality of this tool does not permit to study fine points of solutions to (11.11) . 

Solutions to (11.11) enjoy interesting properties, Fukui and Giga, BFGI . have noticed that facets 
persist. By a facet we mean a flat part (i.e. affine) of the solution with zero slope. Zero is 
exactly the point of singularity of function | • | . This is why the problem of facet evolution is 
not only nonlocal but highly anisotropic. Our equation (11.11) is at least formally parabolic of 
the second order. This is why we call the above behavior of solutions the sudden directional 
diffusion. However, even more dramatic effects of singular diffussion can be seen in the fourth 
order problems, see HGGI 

As we have already mentioned some properties of facets were established in HFGH . e.g. their 
finite speed of propagation was calculated. What is missing is the description of the process how 
they merge and how they are created. In [MR2J we studied a problem similar to (11.11) . We worked 
there with a simplification of the flow of a closed curve by the singular mean weighted curvature. 
We have shown existence of so-called almost classical solutions, i.e. there is a finite number of 
time instances when the time derivative does not exist. However the results of [MR2J indicate 
lack of efficiency of the methods used there. This fact is our motivation to rebuilt the theory 
from the very beginning. For this reason we consider here the model system admitting effects of 
sudden directional diffusion. Hoping that our approach will be suitable for more general systems. 

Our approach is as follows. We notice that the implicit time discretization leads to a series of 
Yosida approximations to the operator on the right-hand-side (r.h.s. for short) of (11.11) . We study 
them quite precisely, because we consider variable time steps. As a result we capture the moment 
when two facets merge. We do not perform any further special considerations. We want to see 
how the regularity of original solutions is transported via solvability of the Yosida approximation. 
Due to the one-dimensional character of the problem we are able to obtain a result so good that it 
is of the maximal regularity character, what is rather expected for quasilinear parabolic systems. 
Let us underline that properly understood smoothness is the most important question connected 
to solvability of the original system. We have to modify standard regularity setting in order to 
capture all phenomena appearing in the system we study. As a result of our considerations we 
come to the conclusion that the best smoothness we could expect for a solution u that u(-,t) be 
piecewise linear function, while x 2 is an example of an irregular function. 

Our main goal is monitoring the evolution, as well creation, of the facets and a precise de- 
scription of the regularity of solutions to (11.11) . which we construct here. For this purpose we 
apply methods, which are distinctively different from those in the literature. We develop ideas 
which appeared in our earlier works. The key point is a construction of a proper composition 
of two multivalued operators: the first one is sgn understood as a maximal monotone graph, the 
other one is u x , which is defined only a.e. We leave aside the issue that in general this is a mea- 
sure, not a function. This problem is resolved differently by the authors applying the semigroup 
approach, ED. HACTH . llGGKl . llBCNH etc. We treat u x as a Clarke differential (see (O) and 
the text below this formula). Here, we show that this composition is helpful when: 

(a) we construct solutions, see Theorem l3.ll 

(b) we discuss regularity of solutions, see Theorems 12.11 and [231 

At the moment, however, the usefulness of this approach is limited to one dimension. The ad- 
vantage of our method is also simplicity, the composition is explicitly computable. As an extra 
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result we obtain asymptotics of solutions. The Dirichlet boundary conditions imply that the set 
of possible equilibria consists of monotone functions. Our analysis shows that steady state must 
be reached in finite time. 

On the other hand, there are two sorts of results available up to now to deal with (11.11) : 

1) the method based on the abstract semigroup theory, see e.g. llFGll . HAC1H . HGGKH and HBCNL 
It is very general and elegant, it enables us to study the facet motion, but it does not capture 
all relevant information. The intrinsic difficulty associated with this method is the fact that the 
energy functional corresponding to (11.11) is not coercive, also see below Lemma I2T1 and the proof 
of Theorem 13. II 

2) the method based on the appropriate definition of the viscosity solution HTGOl . However, a 
different kind of problem was studied there. This is an active research field, see [GGR|. 

Our approach is based on the Yosida approximation, defined by as a solution of the resolvent 
problem 

Xu- — (sgn(u x )) = Xv in (a, b), (u - v)\ 9[aM = 0. (1.2) 

There are a couple of points to be made here. Firstly, we will construct u, a solution to (|1.2I) . 
by very simple means, this is done is Section [31 This process resembles looking for a good 
notion of a weak solutions to a PDE. Since we came up with an integral equation we will call its 
solutions mild ones, see formula (13.161) . 



Secondly, (11.21) may be interpreted as an Euler-Lagrange equation for a non-standard varia- 
tional functional. Namely, we set 

T(„.\ _ / Ja \ Du \ if w £ D(J) = {u e BV[a, 6], u(a) = a b , u(b) = a e }, ( , ~. 

JW I +oo if L 2 (a, b) \D(J), UJ > 

where J a \Du\ is the total variation of measure Du. We stress that we consider the space BV 
over a closed interval. Then, (11.21) may be seen as 

v eu + hdj(u), (1.4) 

where dj is the subdifferential of J and h = \. We shall see that the well-established con- 
vex analysis will yield existence of a unique solution to inclusion (|1.4I) . This solution will be 
called variational. Since variational solutions are stronger (we shall see this), thus both solutions 
coincide. 

We note that the Dirichlet problem in the multidimensional case is much more difficult, in 
particular the meaning of the boundary condition is not clear, see HAC21 

Thus, no matter which point of view we adopted, u is given as the action of the nonlinear 
resolvent operator R(X, A) on v, i.e. 



u = R(X, A)(v) = (X + A)-\v) 



where A = — ^sgn J^. However, the notion of a mild solution to (11.21) does not permit us to 



interpret this equation easily. On the other hand, by convex analysis, we can see (11.21) as an 
inclusion (11.41). 
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The definition of the nonlinear resolvent operator leads to a detailed study of J . One of our 
results, see Theorem 15 .21 is a characterization of solutions to (11.21) . The advantage of (11.21) is that 
it permits to monitor closely behavior of facets. It says that the regularity propagates. That is, 
if v is such that v x belongs to the BV space and the number of connected components of the 
properly understood set {x : v x (x) = 0} is finite, then u x has the same property for sufficiently 
large A. 

It is well-known that the nonlinear resolvent leads to Yosida approximation, which is the key 
object in the construction of the nonlinear semigroup in the Komura theory. Namely, we set 

A x u = X(u- R{X,A)(Xu)). (1.5) 

Our observation is that a maximal monotone multivalued operator like sgn taking values in 
[—1,1] may be composed with a multifunction properly generalizing a function of bounded total 
variation. We shall describe here this composition denoted by o, see Section |2j We introduced 
such an operation in [MR2J, see also [MR3 1. We also point to an essential difficulty here, which 
is the problem of composition of two multivalued operators. Even if both of them are maximal 
monotone, the result need not be monotone nor single valued. If the outer of the two operators 
we compose is a subdifferential, then we expect that the result is closely related to the minimal 
section of the subdifferential. 

One of our main results says that A\u defined by (11.51) indeed converges to — J^sgnou^.. 
Moreover, we have an error estimate, see Theorem 13. 1[ formula (13.41) . In this way we justify 
correctness of the new notion. Due to the "explicit" nature of o, we may better describe the 
regularity of solutions to (|1.2I) . 

Once we have constructed the Yosida approximation, we show existence of solution to the 
approximating problem u x = —A\(u x ) on short time intervals, where u x (t ) is given. This is 
done in Lemma 14.11 In fact, the method is close in spirit to the construction of the nonlinear 
semigroup, see HCLL Convergence of the approximate solutions is shown at the end of Section 
4. Here, we use the full power of the Yosida approximation to capture the finite number of time 
instances when the solution u(t) is just right differentiable with respect to time, otherwise the 
derivative exists. The point is that we can control the distance to the original problem (11.11) . so 
that we can monitor the time instances when facets merge. 

Let us tell few words about the approach of proving our result. First, we define a space 
of admissible functions giving regularity of constructed solutions. Furthermore, we state main 
results together with an explanation of the meaning of almost classical solutions. In Section 3 we 
study the Yosida approximation for our system, concentrating on qualitative analysis of solutions. 
Proofs in this part are based on a direct construction which is possible due to the fine properties of 
chosen regularity. Subsequently, we prove the main results concerning existence and regularity. 
Section 5 is devoted to an alternative proof of existence for the resolvent operator based on 
the classical approach via Calculus of Variations. This analysis shows that restrictions taken in 
Sections 3 and 4 are natural and reasonable. Finally, we study the asymptotics of solutions and 
present an example of an explicit solution. We conclude our paper with numerical simulations. 
They are based upon the semidiscretization. Since they present a series of time snapshots, these 
pictures contain only the round-off error. At each time step there is no discretization error. The 
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examples in Section 6 present the typical behavior, for which each solution becomes a monotone 
function in finite time. 

2 The composition o and the main result 

Our main goal is to present a new approach to solvability of systems of type (|l.ll) . We construct 
a solution with this novel technique, and next we compare it with ones obtained in a more stan- 
dard approach. This will clarify why some assumptions, which seem to be artificial, after deeper 
analysis will look completely natural. The total variation flow is a good example for such exper- 
iment, since we know precisely the solution, additionally its simple form allows us to deduce the 
qualitative properties by standard methods of the calculus of variation. 
The first step is to define the basic regularity class of functions. 

Definition 2.1. (cf. []Zl Chapter 5]) We say that a real valued function it, defined over a closed 
interval [a, b], belongs to BV[a, b), provided that 

\\Du\\ = / \Du\ < oo, 

J a 

where \Du\ is the total variation of the measure Du. We recall that 

IM|w[a,6] = \\Du\\ + ||w||i. 

For the sake of definiteness, but without any loss of generality we assume that [a, b] = [0, 1]. 

Additionally, we treat BV functions as multi-valued function. This is easy for functions 
which are derivatives, u x E BV[a, b). This is very useful in the regularity study of solution to 
(11.11) . Indeed, if u and u x belong to BV[a, b], then u is Lipschitz continuous. Hence, and 

exist everywhere and they differ on at most countable set. Thus, we may set 

d x u{s) = {ru~ + (1 - r)it+ : r G [0, 1]}. (2.1) 

Under our assumptions on it, the set d x u(x) is the Clarke differential of u and equality holds in 
(12.11) due to jCg} Section 2, Ex. 1]. If u is convex, then d x u is the well-known subdifferential of 
it. As a result, if u x E BV, then for each x$ E (0, 1), we have 

d x u(x ) = []jm_u x (x), \im + u x (x)} or , 

where [a, b] or = [a, b] for a < b and [a, b] or = [b, a] for b > a. 

However, the description of solutions as functions whose derivatives belong to BV is not 
sufficient. We have required to restrict our attention to its subclass. There is a need to control 
the facets, which we shall explain momentarily. A facet of u, F is a closed, connected piece of 
graph of u with zero slope, i.e. F = F(^,£ + ) = {(x,y) : y = const = u(x ), x E 
which is maximal with respect to inclusion of sets. The interval £ + ] will be called the set of 
parameters or preimage of facet F. 
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Let us recall that zero is the only point, where the absolute value, | • |, the integrand in the 
definition of J, fails to be differentiable. Thus, the special role of the zero slope and facets. 

We shall also distinguish a subclass of facets. We shall say that a facet F(£~ , £ + ) has zero 
curvature, if and only if there is such e > 0, that function u restricted to [£~ — e, £ + + e] is 
monotone. In the case the function under consideration is increasing this means that — e) < 
u{£~) = < u(£ + + e). We shall see that zero curvature facets do not move at all. There 

may be even an infinite number of them. They have no influence on the evolution of the system. 
For that reason we introduce the following objects, capturing the essential phenomena. We shall 
say that a facet F((~, ( + ) of u is an essential facet. It will be denoted by F ess ((~, ( + ), provided 
that there exists e > such that either 

u is decreasing on (£~ — e,(~) and u(t) > u((~) for t G (£~ — e, £~) and u is 
increasing on (( + , ( + + e) and u(t) > u(( + ) for t G (C + , C + + e ) (then we call such 
a facet convex); moreover we set 

sgnre [c - iC+] = 1 (2.2) 

or 

u is increasing on ((~ — e, £~) and u(t) < u((~) for t G ((~ — e, C~) an d u is 
decreasing on (( + , ( + + e) and u(t) < u(( + ) for t £ (C + , C + + e ) (then we call such 
facet concave); moreover, we set 

sgnK[ C - c +] = -1. (2.3) 

It may happen that £~ = ( + =: (, then we shall call () a degenerate essential facet. In this 
case u has a strict local minimum or a strict maximum at point (. 

We will call sgn the transition number of facet F((~, ( + ). For the sake of consistency 

we set the transition number sgn to zero for a zero curvature facet F((~, ( + ). 

The union of parameter sets of all essential facets is denoted by E ess (w) and K ess (w) is the 
number of essential facets, including degenerate facets. 

Definition 2.2. Let us suppose that w = d x u e BV[0, 1], where u is absolutely continuous and 
is the Clarke differential of u. We define E(w) = {x e [0, 1] : G We say that u; as 

above is J -regular or shorter iu G J-R iff the set E ess (w) C E(w) consists of a finite number of 
components, i.e. 

E ess (w) = [a u 6i] U . . . U [flff esj( „), &x ess ( w )] where Oj < 6j (2.4) 

and each interval [a 1 , b l ] is an argument set of an essential (nondegenerate or degenerate) facet 
F(a\ b % ). In particular, components of E(w) \ E ess (w) consists only of arguments of zero curva- 
ture facets of u. 

Our definition in particular excludes functions with fast oscillations like x 2 sin ^. We distin- 
guished above a subset of BV functions. 
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Since degenerate facets will be treated as pathology, for given w G J-R, we define 

L{w) = min{6 — a : [a, b] is a connected component of E ess (w)}. (2.5) 

Note that L(w) = iff there exists a degenerate facet of u. 

The name J-regular refers to the regularity of the integrand in the functional J , which has 
singular point at p = 0. J-regularity of w = du x means that function u can be split into finite 
number of subdomains where it is monotone. 

We also define the following quantity, 

||H|j-R[ ,i] = IMIw[o,i] +K ess (w), (2.6) 

where K ess {w) is the number of connected parts of E ess (w), however, this is not any norm in this 
space. 

We start with the definition of a useful class of admissible functions. 

Definition 2.3. We shall say that a function a is admissible, for short a G AF[0, 1], iff a : 
[0,1] -^E, 

a = d x a with a G J-R and a(0) = a&, a(l) = a e . (2.7) 
Here, d x a denotes the set-valued Clarke differential of a. 

We note that the above definition restricts the behavior of admissible function at the boundary 
of the domain. Namely, if a G AF, then a is monotone on an interval [0, x ) for some x G (0, 1) 
and either 

a(x ) > a(0) or a(x ) < a(0). 
By the same token, a is monotone on an interval (x , 1] for some x G (0, 1) and either 

a(xo) > a(l) or a(xo) < a(l). 

Thus, the Dirichlet boundary condition makes immobile any facet touching the boundary. 
Hence, such facets behave as if they had zero curvature. 

A composition of multivalued operators requires proper preparations. Due to the needs of 
our paper, we restrict ourselves to a definition of 

sgno a 

for a suitable class of multivalued operators a. Of course, it is most important to define this 
composition in the interior of the domain we work with. See also [MR2J, HMR31I . 

Definition 2.4. Let us suppose that a is admissible and d x (3 = a G J-R[0, 1]. The definition of 
sgn o a is pointwise. Let us first consider x G [0, 1] \ 5 ess (a). Then, there exists an interval (a, b) 
containing x and such that either (5 is increasing on (a, b) or decreasing. In the first case we set 

sgnoa(ar) = l; (2.8) 
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if (3 is decreasing on (a, b), then we set 

sgnoa(x) = -1. (2.9) 

We note that the set [0, 1] \ E ess (a) is a finite sum of open intervals, on each of them function j3 
is monotone. Furthermore, the end points of [0, 1] can not belong to S ess (a). 

Now, let us consider x G E ess (a), then there is [p,q] a connected component of E ess (a) 
containing x. If F(p, q) is a convex facet of /3, then we set, 

2 2p 

sgnoafx) = x 1 for x G [p, q\. (2.10) 

q — p q — p ' 

If F(p, q) is a concave facet of a, then we set, 

2 2p 

sgnoa(x) = x -\ hi for x£\p,q]. (2.11) 

q — v q — v 

We have already mentioned that the Dirichlet boundary condition does not permit any motion 
of the facet touching the boundary. Thus, effectively, they behave like zero-curvature facets. Part 
2 of Definition 12 .41 takes this into account. 

Now, we are in a position to state main results being also a justification of the notion of almost 
classical solutions to our system. 

Theorem 2.1. Let u Q G AF[0, 1], L(u x ) > with Uq(0) = and Uq{1) = a e , then the system 
di.il) admits unique solution in the sense specified by A3.161) and such that 

u x G Loo(0, T; J-R[0, 1]). (2.12) 

Moreover, u is an almost classical solution, i.e. it fulfills di.il) in the following sense 

u t - d x sgaou x = in [0,l]x(0,T), 

u(0,t) = a b , u(l,t) = a e for te[0,T), (2.13) 

u\ t =o = u on [0,1], 

where the time derivative in t\2.13\) exists for all time instances, except for at most a finite number 
of exceptions, the x derivative exists for at most a finite number of exceptions. Additionally, 
u(-,t) G AF[0, 1] fort G [0,T]. 

We study a second order parabolic equation with the goal of establishing existence of almost 
classical solutions. This is why we do not consider general data in L 2 , but those which are more 
natural for this problem, where the jumps in u x and their number matter most. This is why we 
look for u, which not only belongs to BV, i.e. u{-, t) G BV, but also u{-, t) G AF. In addition, 
the necessity of introducing essential facets will be explained. 

An improvement of the above result, showing a regularization effects, is the following 

Theorem 2.2. Let uq be as in previous Theorem above, but L(uo, x ) = 0. Then, there exists 
a unique mild solution to di.il) . which is almost classical and it fulfills A2.13\l . Furthermore, 
L{u x {t)) > Ofort > 0. 
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The second theorem shows that the class of functions with non-degenerate facets is typical, 
and each initially degenerate essential facet momentarily evolve into an nontrivial interval. Fur- 
thermore, creation of such a singularity is impossible. In order to explain this phenomena let us 
analyze the following very important example related to analysis of nonlinear elliptic operator 
defined by subdifferential of (11.31) . 

We first recall the basic definition. We say that w G dj(u) iff w G L 2 (a,b) and for all 
h G L 2 (a, b) the inequality holds, 

J(u + h)-J(u)>(w,h) 2 . (2.14) 

Here (/, g) 2 stands for the regular inner product in L 2 (a, b). We also say that v G D(dJ), i.e. v 
belongs to the domain of dj iff dj(v) ^ 0. 

We state here our fundamental example. We recall (11.31) and for the sake of convenience we 
set (a, b) = (—1, 1). Then we make the following observation. 

Lemma 2.1. Function \x 2 does not belong to D(dJ). 

Proof. If \x 2 G D(dJ), then there existed w G L 2 (-l, 1) such that for all </> G Cg°(-1, 1) 
and t G E 

/ (\x + t(j) x \ — \x\)dx > t / w(f)dx. (2.15) 
J(-i,i) J 

We restrict ourselves to <fi such that 

(f) G C^(-8, 5) and supp0 z [-5, -5/2} U [5/2, 5}. 

Additionally, 

<f) x {t) > for t G (-5, -5/2),</) x (t) < for t G (5/2, -5) 

and 

cj)(-5) = (f)(5) = 0, <p(t) = 1 for t G (-5/2, 5/2). 
Next, let us observe that 

\x + t(f) x (x)\ — \x\ = t(f) x (x)sgnx for 1^(^)1 < S/2; (2.16) 

we keep in mind that 4> x (t) = for t G (—5/2, 5/2). 
Thus, for such and £ the r.h.s. of (12.151) equals 

/ (\x + t(p x (x)\ - \x\)dx = / t(p x ■ (-l)dx + t(f) x - (l)dx = -2t. (2.17) 

J (-6/2,8/2) J (-8,-6/2) J(S/2,8) 

Hence, we get 

—2t>t / wcpdx, 

J(-8,8) 
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what implies for t > 



2 < 



(-5,8) 



wcj)dx < 



(-8,5) 



wldx — > 



(2.18) 




The full description of the domain of the subdifferential dj of (11.31) is beyond the scope of 
this paper. There is a description of D(J) for the multidimensional version of the problem we 
consider, see e.g. HAC2L It is based on Anzellotti's formula for integration by parts HAzl . How- 
ever, a direct characterization of this set for the one-dimensional problem seems to be missing 
even though this functional has been studied in the literature. 

At the end we mention a result describing the asymptotics of solutions, proved in the last 
section. 

Theorem 2.3. There is finite t ext > such that the solution u reaches a steady state at t ext , i.e. 
u(t) = u(t ext ) for t > t ext . Moreover, we have an explicit estimate for t ext in terms ofu , see 



The above result shows that the limit of any solution, as time goes to infinity, is always 
a monotone function, and this will be proved and illustrated in Section 6. There we present 
numerical simulations based on the analysis of system (|l.ll) . It is interesting to note that in 
comparison with [FOP] who deals with the multidimensional case, our computations do not 
contain any discretization error. A rich possibility of stationary states are allowed thanks to 
Dirichlet boundary conditions. Note that such picture is impossible for Neumann boundary 
constraints, for which there are only trivial/constant equilibria. 

3 Yosida approximation 

The central object for our considerations is the Yosida approximation to — c^sgnc^. First, we 
introduce an auxiliary notion of a nonlinear resolvent operator to the following problem, 



where v is a given element of L 2 (0, 1). 

Definition 3.1. An operator assigning to v E J-R a unique solution, u G J-R, to (13.11) will be 
called the resolvent of A = — c^sgn d x and we denote it by 



Am — — sgn (u x ) = \v on [0, 1] 



u = v at d[0, 1], 



(3.1) 



u = R(X,A)v. 



Now, we may introduce the Yosida approximation to A. 
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Definition 3.2. Let us assume that A = —d x sgnd x is as above and A > 0. An operator A\ : 
J-R — > J-R given by 

A\u = X(u- R(X,A)(Xu)) 
is called the Yosida approximation of A. 

Since the notion of Yosida approximation seems well-understood, we will use it to explain 
the meaning of A. For this purpose we will fix w E J-R and A > 0. We set u x := R(X, A)w. We 
will look more closely at A\(u x ). 

Theorem 3.1. Let us assume that w E AF[0, 1], i.e. w x E J-R, then there exists a unique 
solution to 

Xu + A(u) = Xw in (0,1), u(0) = w(0), u(l) = w(l), (3.2) 
denoted by u x , fulfilling 

\\ u x\\bv[o,i] < 1 1 w x 1 1 BV [0,1]- (3-3) 
Moreover, there is A > such that 

K ess (u x ) = K ess (w x ) for A > A with ||u*||j-r < H^Hj-r. 

Furthermore, if L(w x ) = d > 0, equation Q.2D can be restated as follows 

Xu x - d x sgnou x = Xw + V(X,x), (3.4) 

where V(X, x) — > in L q for all q < oo as A — > oo. In addition 

A x (u x ) — > —d x sgnow x in L q (0,l) with q < oo. 

Proof. We would like to present an independent proof of existence of solutions to system 
(13.21) . which is based on simple tools, without any explicit reference to calculus of variations. For 
this purpose, we restrict ourselves to w E AF and for sufficiently large A. A simple construction 
of u x for a given w based upon Lemma IBTTl is presented below. 

Our assumptions give us 

K eaa (w x ) 

Zess(w x ) = |J [<,&:] (3.5) 
i=l 

with a{ < b\. Moreover, a\ > and bf eaM < 1. 

Below, we present a construction of u x . Namely, we consider system (13.21) in a neighbor- 
hood of preimage of an essential facet [a*, b\} of w (it may be degenerate) and we prescribe the 
evolution of this facet. If A is sufficiently large, then we keep the number K ess constant. 

Lemma 3.1. Let us suppose that w satisfies the assumptions of Theorem I3.il Then, for suffi- 
ciently large X, and for each i — 1, . . . , K ess {w x ) there exist monotone functions 

X i — y di(X) and X !->■ &«(A), 



11 



which are solutions to the following problem, 

(V(\)-a l (\))w(a\\))= / w + 2-sgnK [al)H] , iy(6*(A)) = ttf(a*(A)). (3.6) 

TTiese solutions are defined locally, i.e. in a neighborhood of [a* , 6* ]. 

We recall that, the transition numbers sgn/cui^n were defined in H2.2\) . ( 12.31) . Additionally, 
we require 

a x (A) > 0, ^-(^)(A) < 1 and 6* (A) < a m (A) /or i = 1, . . . , i^ ess (w K ) - 1. (3.7) 

However, if A is the greatest lower bound of A as above, then one of the three possibilities 
occurs, 

a 1 (A )=0 or b KM (\ ) = l or a*(A ) = b i+l (X ). (3.8) 
is worthwhile to underline that the lemma holds if L(w x ) = 0, too. 

Proof. Let fix % in {1, ... , K ess (w x )}. Problem (13.61) comes from integration of equation 



(13.21) over a neighborhood of facet [a\ , £>*]. For r G R in a neighborhood of zero and such that 

rsgn/«[ a i i6 i] > 0, we set 

a\r) = min(w| [6 i-i )4] )- 1 (w(<) + r), ^(r) = max^^i^- 1 ^!) + r). (3.9) 

This definition is correct, because functions wLi-i tt< -, and u>Lj a i+i, are monotone. If these 
functions are strictly monotone, then w~ 1 (w(bl) + r) is strictly monotone too, so the min/max 
are redundant. However, if there exists {a} ^ [a, (3} C E(iu) and C [6* _1 ,a*] (resp. 

[a, /3] C [6* , a* +1 ], then (iwLj-i a i]) _1 (resp. (iwLi 4+i]) _1 ) is a maximal monotone graph and 
min/max makes a l (-) (resp. &*(•)) single valued and discontinuous. However, the function 

fb'(r) 

r ^ {b\r) -a l {T))w{a\r)) - w{s) ds =: F^r), i = 1, . . . , K ess (w x ), 

is continuous. Indeed, if r is point, where a % and b % are continuous, then this statement is 
clear. Let us suppose that at r function a 1 has a jump (the argument for b % is the same). Then, 

[a z (r ), 0\ C E(w x ), where a*(r ) < f3 and for any x G [a l (r ), /3] we have 

p¥(T ) _ f¥{r) 

(b\r ) - a i (r ))tw(a i (r )) - / u>(s) = (6*(r ) - s)tw(a;) - / iw(s) rfs. (3.10) 

«/ a 1 (to ) •/ x 

This is so, because we notice that w restricted to [a l (r ), f3] is constant and equal to w(a\) + r . 
Moreover, 

/ w(s) ds = / w(s)ds + / iu(s) (is = (x — a 1 {To)){w(a\ ) + r ) + / w(s)ds. 



'a' (to) Ja l (To) 
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Hence, our claim follows, i.e. continuity of Fi, i = 1, . . . , K ess {w). Indeed , let us suppose 
that r n converges from one side to r (the side, left or right, depends upon sgnK[ alji) ,]) so that 
limbec d l (r n ) = 7, where 7 = a l (r ) or 7 = (3. Then, due to (13.101) we deduce continuity of Fi. 

Subsequently, if we take A sufficiently large, then f sgn «[ a j,6i] is in the range of Fi, i.e. there 
exists Ti = Ti(X) such that F,(r(A)) = |sgn k^^]. If we further make A larger, then we can 
make sure that for each i — 1, . . . , K ess (w x ) we have 

p- 1 (r,(A))<a J (r J (A)) and 6*(r<(A)) < a i+1 ( n (X)). 

Thus, we set 

a l (A) :=a*(r 4 (A)), 6*(A) := ^(A)). 

Let us define A to be the inf of A's constructed above. 
We see that for A one of the inequalities 

^(Ao) > 0, 6 l (Ao) < a m (A ), z = 1, . . . , 4 s ( Wl ) - 1, 6^^)(A ) < 1. 
become equality. □ 

This lemma permits us to define the function u for A > A , 



(w for x G [0,1] \Uf=T ( ^ ) [^(A),& l (A)] n ln 

1 w(a*) for a; e [a* (A), 6*(A)] ^ J 



We immediately notice that K ess (u^.) = K ess (w x ) and S ess (w A ) = U^T^' > [ al > ^1' provided 
that A > A . 

Let us analyze what happens at A = A . We know that one of the three possibilities in 
(13.81) occurs. We notice that if a 1 (A ) = or b Keaa( ~ Wx \\ ) = 1, then a facet of u x touches the 
boundary. Subsequently this facet becomes a zero curvature facet, for it is immobile. This is a 
simple consequence of Dirichlet boundary conditions which do not admit any evolution of facets 
touching the boundary. 

Let us look at the case b z (Xo) = a l+1 (Ao) for an index i. Thus, we obtain the phenomenon of 
facet merging. In both cases the structure of the set E ess (u x ) will be different from E ess (w x ). As 
a result, we have 

K ess (u x x ) < K ess (w x ). (3.12) 

It is worth stressing that at the moment A = A more than two facets may merge, so we can 
not control the decrease of number K. In this case we have to slightly modify (13.111) . since 
the structure of E ess (u x ) is different from E ess (w x ). It is sufficient to notice that the number of 
elements in the decomposition (13.51) has decreased. 

It is clear that for A > A , we have 

K ess (u x x ) < K ess {w x ) (3.13) 
and by the construction, (13.1 II) it is also obvious that (see Definition ^. II) 

\\Du x x \\ < \\Dw x \\. (3.14) 
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Note that the boundary conditions are given, so (13.141) controls the whole norm of 



u 



Once we constructed a solution u x by (13.1 II) . we shall discuss the question: in what sense 



does it satisfy equation (11.21) . One hint is given in the process of construction a*(A) and b l (X). 
This is closely related to ideas in IMRTTl. If we stick with differential inclusions, then formula 

1 d 

-sgnw^. 3 0, (3.15) 



u 



w 



A dx 



leads to difficulties, because we did not provide any definition of the last term on the left-hand- 
side (l.h.s. for short). 

Here comes our meaning of a mild solution: for each x E [0, 1], the following inclusion must 
hold 

1 



w)dx' 



A 



sgnu x 



3 0. 



(3.16) 



We shall keep in mind that at x = 0, we have u = w (for the sake of simplicity of notation we 
shall suppress the superscript A, when this does not lead into confusion). 

In order to show that u fulfills (13.161) . we will examine a neighborhood of the first component 
of E ess (u x ), i.e. [a 1 , b l \. We take x E [0, a 1 ), then u = w on [0, x]. Thus, it is enough to check 
whether ^(sgn u x (0) — sgn u x (x)) 3 0. We notice that on [0, x] C [0, a\) function u is monotone. 
As a result sgn 1^(0) and <Agnu x (x) may equal 1 or [—1, 1], provided that u is increasing. If on 
the other hand, u is decreasing on [0, x\, then sgn 1^(0) and sgnu x (x) are equal to —1 or [—1, 1]. 
If any of these possibilities occurs, then (13.161) is fulfilled. 

We shall continue, after assuming for the sake of definiteness that facet F(a 1 , b 1 ) is convex. 
The argument for a concave facet is analogous. 

Let us consider x E [a 1 , b 1 ]. We interpret sgnt as a multivalued function such that sgnO = 
[—1,1]. Then, we have for x E [a 1 , b 1 ] 



f x 1 1 

J (u - w)dx - -[-1, 1] + -sgnw^l^o 3 0. 

Since we assumed that the facet F(a 1 , b 1 ) is convex, from (13.61) we find that 

< / («- w)dx' < -. 
Jo A 

By the assumption we know that sgn u a 



[u 



x \x'=q 3 -1- Hence, 
w)dx'-\ G^-1,1]. 



(3.17) 



(3.18) 



(3.19) 



This shows (13.161) again. In case F(a 1 , b 1 ) is concave, the argument is analogous. 
Let us now consider x E (b 1 , a 2 }, then we have 



{u — w)dx' — — sgnu x 



{u — w)dx' — — sgnit a 



+ 



(u — w)dx' 



-^sgnw x 



h + h + h- 



6 1 



(u — w)dx' — — sgnw 3 



(3.20) 



b 1 
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Here, we do have the freedom of choosing sgnu x at x = b 1 . Namely we set sgnM 2 .(6 1 ) = — 1. 
We also know that sgnw^a 1 ) = 1 . 

We recall that by the very construction of a 1 and b 1 , we have I2 = 0. Subsequently, we notice 
that the argument performed for x E [0, a 1 ) applies also to x E (b 1 , a 2 ], Thus, 

h + h + h = -^(l-sgn M;E (0)) + 0-^(-l + sgnM x .(x)) 
= -(-sgnu^O) + sgnu x (x)) 3 0, 

i.e. (13.161) holds again. 

Repeating the above procedure for each subsequent facet, we prove that u given by (13.131) 
fulfills (13.161) . The case x E [b Kes3 , 1] is handled in the same way. Thus, we proved the first part 
of Theorem 13 . 1 1 concerning existence. 

We shall look more closely at the solutions when A = Ao- We have then two basic possibili- 
ties: 

(1) The first facet F(a 1 , b 1 ) or the last one F(a k , b k ) touches the boundary, i.e. a 1 = or resp. 
b k = 1. If this happens, then F(0, b 1 ), resp. F(a k , 1), has zero curvature. 

(2) Two or more facets merge, i.e. there are i, r > such that 

lim 6 i_1 (A) = 6 i_1 (A ) < a^Ao) = lim aHX) 

A^Ao A->Ao 

and 

lim b i+j (X) = lim a i+1+j (X), j = 0, 1, . . . , r - 1, 

A— s>Ao A->Ao 

and 

lim b i+r -\X) < lim a i+r (X). 

A->Ao A— s>Ao 

We adopt the convention that b° = and a k+1 = 1. 

When this happens, we have two further sub-options: 

(i) an odd number of facets merge, then F(a l (Xo), b l+r (Xo)) has zero curvature; 

(ii) an even number of facets merge, then [a l (A ), b l+r (X )} 

Of course, it may happen that simultaneously a number of events of type (2i) or (2ii) occurs. 

First let us observe that u = w away from the set {u x = 0}, so we conclude E(w x ) C E(u x ). 
More precisely, the equality holds on a larger set. Namely, if F(a\ b l ) is a zero curvature facet 
and A > Ao, then the very construction of a*(A), b l (X) implies that u = w on [a\ b % ]. If u(a) = 
u(b) = w(a) = w(b), so there must be a point c E (a, b) such that E w x . Thus, we obtain for 
any A > 

K ess {u x ) < K ess {w x ). 

Let L{w x ) — d > 0, then we consider 

Xu + A x (u) = Xw forA>A , (3.21) 
where we suppressed the superscript A over u. 
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As we have already seen taking large A, i.e. A > A , excludes the possibility of facet merging 
or hitting the boundary, thus K ess (w x ) = K ess {u x ). Let us emphasize that K ess (u x ) may decrease 
only a finite number of times. 

Let us suppose that [a*, b*] is a connected component of S ess (u x ), i.e. a* = a*°(A), b* = 
b' l0 (X) for an index i . Without loss of generality, we may assume that this facet is convex. So, 
integrating (13.211) over [a*, b*], we find 

b* rb* 

Am - / Xw = 2. (3.22) 



a" 



First, we want to find an answer to the following question. What we can say about the 
behavior of the following quantity J" + (Xu — Xw), where [a, b] is a connected component 
of E ess (w x ) contained in [a*, b*]. In fact we assume, that a = a l °,b = b l ° . 

Since d = L(w x ) is fixed and positive we find from (13.221) that 

b* rb 

X(u — w) > / X(u — w) > dX(u — tu)|[ a ,&], 

a* J a 

Because u — w is monotone on [a, b). As a result, 

X(u-w)\ [aM <-. (3.23) 

Then we conclude that 

*b* 

X(u -w)< (b* - b)X[w(b*) - w(b)). 

At the same time (|3.23l) yields, w(b*) — w(b) < On the other hand, w is monotone on set 
(b, b*). Hence (T3T231 implies that 

b* -b = b io (X) -b ia <^- 1 (_) J (3.24) 

dX 

where W^ 1 ^) is a strictly monotone (possibly multivalued) function, equal w^ 1 (restricted to an 
interval of monotonicity) plus a constant such that lim t ^ + = 0. Eventually, we get 

b * 2 2 

X(u-w) < )- ^0 as A^oo. (3.25) 

dX d 

Since the analysis for (a*, a) is the same, hence (13.241) and (13.251) imply that 

-b* 

\( u - w ) = 2 + V(X) with V(X) ->• as A -> oo. 

Note that V(X) depends only on w, so in Section|4]we will study the approximation error V(X) 
and we will show uniform bounds, provided that L(w x ) > d > 0. 
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Integrating (13.211) yields 



b* pb* 
X(u-w) = -A x (u) = 2, (3.26) 



a-' 



but the pointwise information from the equation yields 

X(u - w)\[ ajb ] = -A x (u) = const. (3.27) 
Thus, taking into account (13.261) and (13.271) . we get 

— Aa(w)|[ ,6] — > 2/(6 — a) a.e. as A — > oo. 

Here, we used that a* = a 10 (A) —> a 10 , b* = b l0 (X) — >• b l ° as A goes to infinity. But a, b depends 
only on w, additionally we shall keep in mind that (13.231) via (13.211) implies that < 2/d 

on whole [0,1]. 

Clearly, by Definition 12 .41 

2 

^sgno^ = — — for x G [a*, b*}. 

[b* — a* ) 

Hence, we have proved that 

A x (u) = -d x sgnou x + V{X,x), (3.28) 

where V(X) = j\ V(X, x) dx and V(A, x) — > in at least L\{I). Here, we should note clearly 
that all depend on A, since a* = a l °(X), b* = b l °(X). We see that we have already proved that 
\V(X,x)\ < |, and/x({suppy(A,.)}) -»■ which gives a relatively strong convergence. Note 
that in (|3.28l) we are not able to obtain "pure" discontinuity in the composition 6, since we work 
with solutions only, hence sgn ou^. must be piecewise linear. 

Next question is: whether c^sgn ou^. — > — <9 x .sgn ow x and in which space? 

Let us observe that (see Definition 12 .11) 

H-Dw^ll < Ij-DiUarll and u^. — )■ w x in measure on /. (3.29) 

It follows that 

sgno Ux -\s( Wx ) ->■ sgno^lg^ uniformly. 

We remember that sgn ou^ and sgn ow x are piecewise linear functions and the set S(w x ) is inde- 
pendent from A, but the case L(w x ) = d > implies that 

A(u x ) -> -—sgnow x inLy(0, 1) qe[l,oo). (3.30) 
Theorem l3.1l is proved. □ 



In particular, as a result of our analysis, we get that the constructed solution to (13.21) is varia- 
tional. 



17 



Lemma 3.2. Function u x given by Theorem \3.1\ is a variational solution to ( 13.21) . i.e. u x fulfills 
(Aw\ 0) + (a(x), 0') = (\w, 0) for each G C °°(0, 1) (3.31) 
and a(x) G sgn o u x (x), where here o denotes the standard composition. 
Proof. From the inclusion (13.161) . we are able to find such a that 

/ (u-w)- \<t(x) + \a(0) = 0. (3.32) 
Jo A A 

Then, testing it by 0' with G C^°(0, 1), we get (13.311) . In particular, we already have shown 
that \R(\, A)X is a monotone operator in L 2 . □ 

4 The construction of the flow 

A key point of our construction of solution is an approximation of the original problem based 
on the Yosida approximation. Here, we meet techniques characteristic for the homogeneous 
Boltzmann equation HdBllMll . For given A, t and A\ defined in (11.51) . we introduce the following 
equation for u x , 

/•to+t 

u x (t + t ) = u x (t )- A x (u x )ds, w A (0,t + t) = a b , u x (l,t + t) = a e for t G (0,T). 

J to 

(4.1) 

We stress that its solvability, established below, does note require that L(u x (to)) > 0. 

Lemma 4.1. Let us suppose that u x (t ) G J-R(J), where I = [0, 1], then there exists a unique 
solution u x to ( 14.71) on the time interval (t , t Q + and 

u x eC(to,t +—;L 2 (I)). 

Moreover, 

sup \\u\t + t)\\j- R < ||u A (t )||j- R . (4.2) 

*e(o,3x) 

Proof. We will first show the bounds. Let us suppose that u x is a solution to (14.11) . then 
Definition [32] and the observation j- t [e xt u x ] = -e xt A x (u x ) + Xe xt u x imply that, 

/■to+t 

u x (to + t)= e- xt u x (t ) + / e- A(t0+ ^ s) A J R(A, A)Am a (s)c/s. (4.3) 

J to 

In order to obtain the estimate in BV, we apply Theorem l3.1[ inequality (13.31) . getting 

su P ||m a || w < e- xt \\u x {t )\\ BV + sup\\R{\,A)\u x (t)\\ BV [ Ae" As rfs 
t t Jo 

< e- Ai ||^(t )IUv + sup ^||An A (t)|| w (l - e~ A *). 
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So we get 

sup IIm^Ubv < \W^{to)hv- (4.4) 
t 

In order to prove existence, we fix A (we will omit the index A in the considerations below) 
and we define a map $ : (7(0, T; L 2 (I)) -> (7(0, T; L 2 (-0) such that v = $(w), where 

v(t) = e-**v + f e x{t ~ s) XR(X,A)Xwds. (4.5) 
Jo 

We notice that due to E((XR(X, A)Xw) x ) D E(w x ) we obtain E(v 0)X ) C S(w x (t)) for t G (0, T), 
provided that w\ t =t = v . Combining this observation with w\ t =t Q = vq again yields, 

E(v 0iX )cE(v x (t)) for te(0,T). (4.6) 

We see that a fixed point of the above map yields a solution to (|4.1I) after a shift of time. For 
the purpose of proving existence of a fixed point of $, we will check that $ is a contraction. 
We notice that if w, w G (7(0, T; L 2 (I)), then monotonicity of i?(A, A)X (thanks to Lemma 13721) 
implies that 

\\R(X, A)Xw - R(X, A)Xw\\ L2 < \\w - w\\ L . 2 . 

Hence, 

||$H - 11^(0,1^(7)) < £ Ae^- s )d S || J R(A, A)Au; - R(X, A)Xw\\ Lao{0>T . M{I)) 

< (1 - e" TA )||w - U)||L oo (0,T i L a (J)), 

i.e. $ is a contraction provided that < T < Now, Banach fixed point theorem implies 
immediately existence of u x , a unique solution to (|4.1I) in (7(0, T; L 2 (I)). 

An aspect is that the solution to (14.31) can be recovered as a limit of the following iterative 
process 

v k+i = (47) 

We have to show that the fixed point belongs to a better space. For this purpose we use 
estimate (1441) . which shows also that if = M, then ||^||sy < M for all k G N. 

Moreover, convergence in L 2 (I) implies convergence in Li(I) and lower semicontinuity of the 
total variation measure (see [Z, Theorem 5.2.1.]) yields u x G L^O, T; BV(I)). 

Finally we show that 

K ess (u x (t + -L)) < K ess (u(t )). (4.8) 
For this purpose it is enough to prove that 

u x (t + t) =u x (t ) on I\E(u x (t +*)) foralU < — , 
but Theorem [3TT1 implies 

i?(A, A)Xu x = Xu x on I\E(R(X,A)Xu x ), 
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namely A x (u x ) = at / \E(R(X, A)Xu x ). Additionally (1461) yields that E(u x (t )) c S(w A (t + 
what finishes the proof of (14.81) . 

Thus, the definition of the solution to (|4.1I) as the limit of the sequence (14.71) together with 
(14.81) imply (|4.2I) . The Lemma is proved. □ 

Lemma 4.2. Let us consider u x (-) given by Lemma \4.1\ IfL{u x (to)) = 0, then L(u x (to + ^)) > 
0. 

Proof. Let us assume a contrary, then there exists a degenerate facet F[a\ b % ] with a* = b l 
such that all functions u x (t + 1) are convex in a neighborhood (p, q) of point a* and they all have 
a minimum only in point a 1 . If functions u x (t + 1) are concave, then the argument is analogous. 
Let us then integrate (14.11) over (a', b') such that a' 1 G (a', 6') C (p, g), 

6' /-fe' rto+t rb' 

u x {t + t)= u\t )- / A\(u x )ds. 

' Ja' Jt Ja> 

But 

rb' rb' 

A x (u x )= \(u x - R(\,A)\u x ) = -2, 
because m a is convex on (a', 6'). Hence, we find 



b' r b' 

A/j- , jA / „.A/ 



u x (t + t) = / U A (t ) + 2t. 

Ja' 

But if our assumption that a 1 = b l were true, then we would be allowed to pass to the limits, 
a' — > a 1 and b' — > a 1 concluding that = + 2t, which is impossible for positive t. Thus, we 
showed that u x (t + ^) does not admit degenerate facets. □ 

After these preparations, we finish the proofs of Theorems 12. l\ and \2. 21 We shall construct an 
approximation of solution on a fixed time interval, say [0, 1]. Let us assume that 

U x : [0, 1] x / R 

is given as follows 

k k + 1 

U x = u x k for t G [— , — — ) and < k < 3A, 
oA oA 

where functions {u x } are given by the following relations 

ft 1 
u x (t) = u - J A x (u x )ds for te (0, - 

M A (t! + t) =M 1 (t 1 ) 



3A J 

ti+t 



jT 1 A A ( M A )rf S for t G (0, 
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u X k+i(tk +t) = u k {t k ) - / A x {u x k+1 )ds for t e (0, — ], 
«ax(*3A-i +t) = wJx_ a (t3A-i) - y A A («^)da for * e (0, — 



and t fc = ^ for < k < 3A. 

||^ A ||l oo (0,T;J-R) < ||Wo||j-R. 

Now, we pass to the limit with A. The estimates imply that \\U x \\l co (o,T;L 2 (i)) — M. Thus, 
we can extract a subsequence such that 

U x ^* U* weakly * in L oo (0, 1; L 2 (I)). 

Moreover, the lower semicontinuity of the total variation measure yields 

\\U x {t)\\ BV < ||w(0)||bv fora.e.t G [0,1]. 

Thus, we should look closer at the limit 

rta+t 

U*(t + t) = U*(t ) - lim / A x (U x (t + t)) ds. 

A ^°° J t 

Let us observe that for a fixed A the function K ess (U x (t)), taking values in N, is decreasing, 
so facet merging may occur just only a finite number of times. 
Let K{uq) = k°, then for a given A we define T x as follows 

K ess (U x (t)) = k° for te[0,T?) and K ess (U x (T x )) < k°. (4.9) 

For a subsequence \imT x =: T\. Indeed T x = T x ' for all sufficiently large A, A' see Lemma l54l 
so we have here T\ > 0. However, we prefer to consider a more general argument valid for more 
complex operators. 

In a similar manner to (14.91) we define a sequence of time instances {T fc }^L 1 . By the defini- 
tions, for any e > there exists A e , such that for A > A e - up to possible subsequence - we can 
split the time interval [0, 1] into following parts 

[0, 1) = [0,Ti - e) U [1\ - e,T 2 + e) U [T 2 + e,T 3 - e)] U ... U [T m + e, 1) 

and 

K ess (u x (t)) = K ess (U*(t)) for te[T k + e, T k+1 - e), 

so {Tfc} is a finite sequence of moments of time at which facets merge. In order to avoid unneces- 
sary problems we restrict ourselves to a suitable subsequence guaranteeing the above properties. 
Now. Theorem 13.11 yields A X (U X ) — > A(U*) = — d x sgn o U* in L q (Q, 1) on time intervals 



(T k + e, T k+ i — e), since by (13.281) we control this convergence uniformly at whole intervals. So 
we get 

U*(t + t) = U*(t )- [ A(U*(s))ds, 

J t 
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because we consider one interval [T k + e, T fc+1 — e). However, crossing T k requires some extra 
care. 

In order to extend the result on the whole interval [0, 1], it is sufficient to prolong the solution 
onto interval [T& — e, T k + e). For this purpose we can use that u x belongs to C(0, 1; L\(I)), see 
Lemma |4~T1 Continuity of of the solution allows us to cross points T k . It follows that 

-fU* exists except points {T k } 

and by the properties of solutions on intervals [T fc , T fc+1 ) we find that the right-hand-side time 
derivative exists everywhere, including points {Tk} 

~jr^U* exists everywhere on [0, 1]. 
Finally, we have shown that U* fulfills 

w u ' = -l s ^ u " (410) 

as an almost classical solution. 

By construction U*(t) E AF. additionally Lemmal4~2lyields L(U*(t)) > Ofort > 0, even as 
L{uQ tX ) = 0. Moreover, the features of almost classical solutions imply that they are variational, 
too. Hence, the monotonicity of sgn implies immediately uniqueness to our problem. Theorems 
12.11 and 12.21 are proved . □ 

Now we want to obtain the same result starting from the classical point of view of the calculus 
of variation in order to explain the chosen regularity. 



5 The variational problem 

In this section we will prove Theorem 13.11 using the tools of the Calculus of Variations. This 
result establishes existence of solutions to (13.11) . i.e. 

Am — 7-sgn (u x ) = Xv in (0, 1), u = v for x = 0, 1 

ax 

for an appropriate v. 

Some parts of the argument, when v G J-R with L{v x ) > are a repetition of results from 
Section [31 However, this repetition is necessary in order to explain that approach from previous 
sections are based on a reasonable class of function, which can be viewed as typical. 

It is clear that first we have to give meaning to this equation. We can easily see that it is 
formally an Euler-Lagrange equation for a functional J h)V defined below. 

Jh,v(u) = hJ(u) + ^ ( (u-v) 2 , 
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where J is introduced in (11.31) . When no ambiguity arises, we shall write J v in place of Jh,v 

We notice that J v is proper and convex. Momentarily, we shall see that it is also lower 
semicontinuous, hence its subdifferential is well defined, see flBrJ, in particular D(dJ) ^ 0. We 
recall that u E D(dJ) if and only if dj(u) ^ 0. It is a well known fact that u is a minimizer of 

J v iff 

hdj(u) + u - v 3 0. (5.1) 

Since hdj{-) + Id is maximal monotone, then for any v E L 2 there exists u E D(dJ) satisfying 
(I5TTT) . see flM 

In this way, we obtain our first interpretation of (13.11) as a differential inclusion. This is not 
very satisfactory as long as we do not have a description of the regularity of the elements of 
D(dJ). We note the basic observation and present its direct proof. 

Lemma 5.1. (a) For any v E L 2 (a, b) functional J v is lower semicontinuous in L 2 . 

(b)Ifv E L 2 (a, b), then there exists u E D(J) C BV(a, b) a unique minimizer of J v . Moreover, 

\\Du\\ = [ \Du\ < \B - A\ + \ [ (v-£) 2 dx, 
Jo 2/i J 

where i is an affine function such that £(a) = A, £(b) = B. 

Proof, (a) Let us suppose that {u n } C L 2 is a sequence converging to u in L 2 . If 

liminf ||-Dw n || = oo, 

n— ¥oo 

then there is nothing to prove. Let us suppose then that sup ngN ||-Dw„|| < K. By the lower 
semicontinuity of the BV seminorm, we infer that u E BV and \\Du\\ < K. The problem is to 
show that the limit u satisfies the boundary conditions. 

If v E BV[a, b], then there is a representative such that \\Dv\\ = V£(v). Moreover, ess sup \v\ 
is finite, see [|Zl Chapter 5]. Thus, there is a representative v satisfying the boundary conditions 
and V£(v) < \\Dv\\ + 4||t>|| 00 . As a result, we select a sequence of representatives u n satisfying 
the boundary conditions and with uniformly bounded variations. Since u n is a sequence of 
bounded functions with commonly bounded total variation we use Helly's theorem to deduce 
existence of subsequence {u nk } which converges to u°° everywhere. Since all functions {u nk } 
satisfy the boundary data, the pointwise limit will satisfy them too. Moreover, due to uniqueness 
of the limit u°° = u a.e. thus we can select a representative belonging to D(J) as desired. 

(b) By definition J v is bounded below. Let us suppose that {u n } is a minimizing sequence in 
L 2 . Of course u n 's belong to BV(a, b) and 



/ \Du n \ + - / (u n - v) 2 dx < K. 
Jo * Jo 



i.e. the sequence {u n } is bounded in the BV norm. Since sets bounded in BV are compact in 
any L p (0, 1), p < oo, see [ABu], we deduce existence of a subsequence {u Uk } converging to u. 
Because of part (a) we infer that u E D(J) and 

|£>u| < liminf I \Du n . \ < K. 
i, k ^°° Jo 
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Combining this with strong convergence of {u nk } in L 2 we come to the conclusion that u is a 
minimizer of J v . 

Uniqueness of a minimizer is a result of strict monotonicity of the operator Id + hdj '. 

Since, u is a minimizer, then J v [t) > J v {u), where £ is an affine function such that £(a) = A, 
£(b) = B. Hence, the desired estimate follows. □ 

We shall establish how much of the smoothness of v is passed to u. Here is our first observa- 
tion. 

Theorem 5.1. IfvE Wl (a, b) , where p E (1, oo) , then u the unique minimizer of 



J v (u)= I h\u x \ + -(u-vy = hj(u) + J -(u-v) 



2 



belongs to Wl and \\u\\i tP < \\v\ 



hp- 



We want to look at the propagation of regularity, so the assumption v x E BV is natural from 
many possible view points. So here is our main result, it will be shown after Theorem 15. 1[ Its 
proof follows from the analysis of the argument leading to Theorem 15. II 

Theorem 5.2. Let us suppose that v e AC[a, b] and u he the corresponding minimizer of J v . 
Then, 

(a) K ess (u) < K ess (v); 

(b) ifv x G BV and K ess (v) is finite, then u x G BV and \\u x \\ B v < H^sIIbv- 

We see from its statement that a type of regularity which propagates is defined by v x E BV 
and a finiteness of the number K ess {v). At this point, we do not claim that this is optimal. 

In order to provide a proof of Theorem 15. 11 we will proceed in several steps. First we shall 
deal with continuous piecewise smooth functions, then we shall show that our claim is true for 
any v which may be approximated in W\ by such functions. We need a simple device to check 
that a function is indeed a minimizer. 

Lemma 5.2. Let us suppose thatv, u E AC[a, b] with v(a) = u(a), v(b) = u(b) and there exists 
o E W\[a, b) and such that o(x) E sgn(u x .(a;)) with sgn understood as a multivalued graph, 
which satisfies the equation 

h—o = u — v (5.2) 

dx 

in the L\ sense. Then, u is a minimizer of J v . 

Proof. Let us take any ip E C^°. Let us calculate 

/>b />b rb -t 

J v (u + if) - J v {u) = h \u x + cp x \-h \u x \+ / -[(u + ip - v) 2 - (u - v) 2 } 



a 



\u x 


+ ¥x\ 


- h 


/ Kl 

/ a 

r\ i 


\u x 


+ (Px\ 


- h 


/ \ u x\ 


\u x 


+ <fx\ 


-h 


r\ i 



2' 

-b 

> h I \u x + (p x \—h I \u x \+ I (u-v)cp 

a J a J a 

// / \u x + <p x \ - h [ \u x \-h [ v-^P- 
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We used (15.21) and the integration by parts. We deal separately with the sets {u x > 0}, {u x < 0} 
and {u x = 0}. We have, 



(J v (u + ip) - Jviu))^ 1 > 

+ 

We used here the fact that a(x) G [— 1, 1] as well. 

Now, we deal with general ip G BV such that u + (p G D(J). We proceed by smooth 
approximation (p n such that tp n converges to ip in L\ and ||-D(^ n || — > \\D<p\\. By what we have 
already shown, we have 

J v (u + (f n ) > J v (u). 

Hence, the inequality is preserved after a passage to the limit. Our claim follows. □ 
We may now start the regularity analysis. 

Lemma 5.3. Let us suppose that v e C[a,b], v(a) = A, v(b) = B, and its derivative exists 
almost everywhere and it is piecewise continuous, its one sided derivatives exist everywhere and 
the sets {v x > 0}, {v x < 0} are open and the number of essential facets of v is finite. Then, for 
any positive h and u a unique minimizer of Jh, v , we have u G Wp with 

\\u\\i, P < IMkp- 

Moreover, there exists a G W^, such that for all x G [a, b] we have cr(x) G sgn (u x (x)) and 
equation (5.2\) is satisfied everywhere except a finite number of points. In addition, 

,, ,, !„ ,, 

F l,oo < 1 + r B oo- 

h 

Proof. We shall proceed by induction. We first show, however, a slightly stronger result if 
v is monotone i.e. the number K ess is zero, and to fix attention we assume that it is increasing. 
Namely, we claim that in this case u = v. We have to show that for any ip such that v+(p G D(J), 
i.e. (p is zero at the ends of [0, 1], we have 

JhA v + V 9 ) > JhA v )- 

Let us notice that 

Jh,v(v + <p) 



We may also set a = 1, since v is increasing. 



/ (| War + Vx\ - U x - 1 • (p x ) 

J{u x >0} 

/ (| Was + <fx\ + u x + 1 • ip x ) + / (\<p x \ - cr ■ tp x ) > 0. 



1 

(h\v x + <p x \ + -tp 2 ) > I h(v x + Lp x ) 



B-A 



hv x = Jh, v (v). 
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The first non trivial case occurs when we have a single essential facet F ess (a, b). The set 
[0, 1] \ [a, b] consists of exactly two components E + (v ) and E~(v). They are such that v\ E +r v \ 
is increasing while v\e-m is decreasing. We stress that the endpoints 0, 1 cannot belong to any 
essential facets. For the sake of fixing attention, we may assume that for all x G [a, b] function 
v has a maximum at x , vm = maxt> (x) = v(x ). We can find £~ G E~(v), £ + G E + (v), i.e. v 
increasing on a] while it is decreasing on [b, and such that 

«(D = Ht) = v com (5.3) 

and v com is the smallest number with this property. In addition, since v is not strictly monotone 
on E + {v) or E~(v), we require that if ( G E + {v) (respectively, £ G E~(v)) is another number 
satisfying (15.31) . then £ < £ + (respectively, £ > In this way £ + , £~ are uniquely defined. 

We want to solve (15.21) . for this purpose we will utilize results of Lemma I3T1 In the present 
case the term — ^ is used in place of A{u). Since we are dealing with a single facet we may be 
more specific about the range of r appearing (13.91) . We notice that for any r G (0, Vm — v com \ 
there exist £~(r) G [£~,a] and£ + (r) G such that 

f + r = v(i + ) + t = v com + r. 

Here, we change the notation and we write £~(r) (respectively, £ + (t)) in place of a(A) (respec- 
tively, 6(A) and h = 1/A. 

In order to solve (15.21) . we have to find simultaneously u and cr(x) G sgnw^., where sgn is 
understood as a maximal monotone graph. We want that u be constant equal to v(£~(h)) on yet 
unspecified [£~(h), £ + (/i)] containing x . On this interval u x will be zero and cr(x) g sgnO will 
be different from zero. Integration of (15.21) over £~(/i), yields an analogue of (13.61) . i.e. 

-2h= (u-v)dx. (5.4) 

J(-(r{h)) 

In Lemma [3TI we established continuity of the mapping 

[0, T: max .) 3 t ^ {v(Ci T )) - v(x)) dx, 

(where T ma x = — v C om)- Moreover, it is strictly decreasing and equal to zero for r = 0. 

Hence, for a fixed h there is at most one r(h) such that (15.41) . If there is such r(h), then for 
the sake of simplicity we shall call ^ ± (t(/i)) by ^(h). Thus, we set 



u[x) 



v(£-(h)) forxe[r(h),e(h)], 
v(x) elsewhere. 



Of course, we have the estimate \\u x \\ p < \\v x \\ p for any p G [1, oo]. 

We have to define a G sgn (u x ). On the set {u x > 0} U {u x < 0}, there is no problem for we 

put 

a(x) = sgn(u x (x)). 
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Before we proceed with the inductive step we introduce a new notation. Let us suppose that 
F(ai, bi), . . ., F(a N , b N ) are all essential facets. Let us look at [0, 1] \ {jf =1 [di, h] consisting of 
open sets (in [0, 1]) (pj, qj), j = 1, . . . , N + 2. Each of the intervals (pj, qj) has the following 
property, either i>|(p-, g -) is increasing, then we write (pj, qj) G E + (v), or v \(p-, q ) is decreasing, 
then we write (pj, qj) G E~(v). We note that the intervals (pj, qj) are maximal sets (with respect 
to set inclusion) with the above property. 

By the very definition, for u as in the statement of this Lemma, we have the following de- 
composition into disjoint sets, 

[0, 1] = E + (u) U BT(u) U ~ ess (u x ) U (E(u x ) \ E ess {u x ). (5.5) 

In general, if u G AF we say that x G E + (u) (resp. x Q G E~(u)), iff x & E ess (u x ) and there is 
(a, (3), a connected component of {u x > 0}, such that there is (1(a), r((3)) containing (a, (3) and 
maximal with respect to set inclusion such that u\ai a ), r [B)) is increasing. In a analogous manner 
we define E~(u). We notice that E + (u) and E~(u) are open and disjoint. We notice that E + (u) 
and E~ (u) are open and disjoint. It is obvious that the decomposition (15.51) is valid for smooth 
u. Moreover, it is not difficult to notice (we will not use it) that if u G AF, the decomposition 
(T531) holds. 

We note that {u x > 0} C E + (v) and {u x < 0} C E~(v) with the possibility of strict 
inclusion. We set a equal to 1 on E + (v) \ {u x > 0} and a equal to —1 on E~(v) \ {u x < 0}. 
Otherwise we define a so that (15.21) holds, e.g. on [£~(h), £, + (h)} we set 

1 f x 

a( x ) = l + - (v(C(h)) - v(x))dx. 

The complement of E + (v) U E~(v) U \^(h), i + (h)} is easy to consider and left to the reader. 
We also mentioned the possibility that 

r£, + {Tmax) 

I / WV(t)) - v{x)) dx\ =: 2h max < 2h. (5.6) 

If this happens we proceed as follows. We find u by the above procedure yielding a minimizer 
of the functional Jh max ,v By Lemma [5~4l we split the minimization problem into two: one 
for Jh max ,v already accomplished and for Jh-h max ,u- Let us notice that the process above for 
h = h max yields u which is monotone. We have already noticed that if u is monotone, then the 
unique minimizer of Jh~h max ,u is u itself. 

Here comes the inductive step. We construct u for v with N + 1 essential facets, denoted 
as above, provided that we know how to deal with v with N essential facets. For each essential 
facet F(cii, bi), i = 1, . . . , iV + 1, we may find intervals i = 1, . . . , iV + 1, constructed 

as above. We may assume that the ordering is such that the sequence of numbers f* \v(x) — 

v(£~) | dx, i = 1, . . . , iV is decreasing. By the process described earlier, for a given positive h, we 
define intervals [^(h),^(h)]. We have two cases to consider: (a) interval Cjv+i(^)] ^ s 

contained in [0, 1] and it does not overlap any of the intervals [^(h),^(h)}, i = 1, . . . , N, i.e. 
£tv+i CO 1S positive, and it is bigger than (h) for all j such that > £~ (h) ; at the same time 
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Cv+iW < 1 an( ^ Cjv+i < £fc (^) f° r ai l k such that < ^(h); (b) the previous condition 
does not hold, i.e. interval [£# +1 (/i)> Cv+iW] * s not contained in [0, 1] or it intersect at least one 
interval [^(h),Ct(h)]. 

The first case presents no problem. The intervals [0, ^ + i{h)), {£,N + i(h), 1] contain no more 
than N essential facets F(a&, Thus, by the inductive assumption we know how to resolve 
any possible overlapping. 

If (b) occurs, then there is j , such that [€j~ (h) , £j~ Q (h)] intersects [^ N+ i(h),^ +l (h)] or 
[£jv+iW> £iv+iW] * s not contained in [0, 1]. The second case is easier, we shall deal with it 
first. It means that there is h < h such that £#+1(^0) — or £^+1(^0) — 1- But then, as we 
know, F(0, £jy +1 (/io)) or F(^ +1 (h ), 1) are not essential facets, thus we consider the minimiza- 
tion of Jh 0!V with minimizer u having N essential facets (of course we have to adjust the integral 
of integration in the functional). If it is so, then by the inductive assumption we are able to re- 
solve any interactions, i.e. intersections of iV essential facets. Then, we solve the minimization 
of Jh-h ,u where the minimizer has no more than N essential facets. 

Thus, inevitably we deal with interactions of facets. Resolving the interactions is easier with 
Lemma l5~4l below, which says that h may be split, if necessary, when ^(hi) = ^(h{), and 
hi < h, while £j < £+. Let us assume that hi is the smallest with this property. We solve 
our problem with v and hi, we find a minimizer of Jh ljV . We may do so, because of lack of 
interactions, we denote its solution by u 1 . Due to the occurrence of interactions the number 
of the essential facets F(a' i , of u 1 is smaller than for v. Thus, we may use the inductive 
assumptions to continue, i.e. to solve our problem with data u 1 and hi — h — hi, in place of h. 
By Lemma [541 solution u 2 is what we need. The proof of the lemma is complete. □ 

Our next Lemma explains that h may be split into smaller steps at will. This permits to 
perform additional analysis at the intermediate steps. 

Lemma 5.4. Let us suppose that v is absolutely continuous and hi, hi > the sets {v x > 0}, 
{v x < 0} are open and they have a finite number of connected components. If u 1 is a minimizer 



of 




while u 2 is a minimizer of 




then u 2 is a minimizer of 




with h — hi + hi. 



Proof. In fact due to our assumptions we have solutions to the equations 



hi— a = u — v 

ax 



hi— a 2 = u 2 — u . 



(5.7) 
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We note that the sequence of implications: u 2 , is different from zero at x, then u\ has a sign there, 
hence v x has a sign too. Moreover, if v x — on an interval (a, 0), then m 2 , are zero (a, /3) 
too. 

We want to show that 

/i— sgn-u?. = -u 2 — i> (5.8) 

dx 

has a solution. Let us add up the two equations above. This yields, 

d (h x x h 2 2 \ 2 

Of course o := ^-a 1 + ^-a 2 G [—1, 1]. If at x we have u x {x) > 0, then v x (x) > 0. Hence, 

/ x hi 1 h 2 2/ \ hi h 2 
<x) = -a (x) + -a (a; ) = _ + _ = l. 

The situation is similar if u 2 x (x) < 0. Let us suppose now that u 2 x (x) = 0, then regardless of the 
sign of ul(x), we know that a{x) G [—1,1] and, by the definition of a, equation (15.81) is satisfied. 
In particular, 

—2h = h / a(x) dx — / (v(£ 2 ) — v(x)) dx = / (u 2 (£, 2 ) — v(x)) dx. 

□ 

The value of this result is that it permits us to split h. We may say that this shows the 
semigroup property. Finally, we show that functions with finite number of essential facets are 
dense in the topology of W 2 X - 

Lemma 5.5. Ifv is smooth with v(a) = A, v(b) = B, then there exist Vk satisfying the assump- 
tion of Lemma 15. 31 Moreover Vk converges weakly tov in W\ and ||vjb||i2 < IMIi,2- 

Proof. The sets E + (v), E~(v) consist of at most countable number of open intervals, 

E ± (v) = {Jl k h (v). 

Subsequently, we suppress the ± superscripts. 

We order the intervals Ik, k G N in so that | If.\ > On U,- =1 Ij, we set v k (x) = v(x). 

On the complement, we define v k to be piecewise linear and continuous. We immediately notice 
that 

ll^lh < K-|| 2 , 

because the linear functions are harmonic. Hence, they minimize the functional J \v x \ 2 with 
Dirichlet data. We have to show that v x converges to v x in L 2 . 

We will show first the pointwise convergence of v k . Let us take any x G [0,1]. If x G 
E + {y) U E~(v), then x G Ij , hence v k (x) = v(x) for k > jo- We suppose now that x is in the 
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complement of E + (v) U E~(v). For the sake of further analysis, we set T k = [0, 1] \ [j i=0 h- 
Each of the sets JF k consists of a finite sum of closed intervals and x e [a k , j3 k ], k e N. By 
construction the sequence a k is increasing, while is decreasing. We shall call by a and /3 
their respective limits. Of course, we have that v k (a k ) = v(a k ) thus this sequence convergence 
to v(ct), while v k (/3 k ) = v(/3 k ) converges to v(/3). We have two case to consider: 1) a < f3, 2) 

a — j3. In the first case we have v k = V ^ k j — v ( ak ) ^ mus t converge to zero. Otherwise, 

Pk - ct k 

we had v x ^ on a subset of (a, (3) of positive measure which is impossible by the definition of 
jys. Hence, v(a) = v(x) = v(/3) for all x e [a, i.e. f fc (x) converges to v(x). 

If a = /3, then our reasoning is similar and by continuity of v we deduce that v(a) = v(x) = 
v(/3). 

Thus, we have shown that v k converges everywhere to v. On the other, hand the bound 
\\v k \\2 < \\ v x\\2 implies that we can select a weakly convergent subsequence. Due to uniqueness 
of the limit it must be v. Since any convergent subsequence of v k converges to u, the whole 
sequence v k converges to v. 

Moreover, due to Sobolev embedding, we deduce that v k converges to v uniformly. □ 

Lemma 5.6. If a sequence v k e W\ converges to v in W\ and u k e W\ is the sequence of 
corresponding minimizers of J v k, then u k converges to u weakly in W\ and strongly in L 2 . 
Moreover, u is a minimizer of ' J v and \\u\\ 1:2 < ||f||i,2- 

Proof. The convergence of u k in L 2 follows from the monotonicity of the subdifferential. 
Indeed, since u k is a minimizer, then 

hdj{u k ) + u k - v k 3 0, 

i.e. there exists ( k e dj(u k ) such that for any test function e L 2 we have, 

MC fc ,0) + (^,0) = (^,0). 

Once we take = u k — ui, we can see that 

MCfc - Cuu k - u') + - u l \\\ = (v k - v l ,u k - u l ). 

Due to monotonicity of dj this implies that \\u k — u l \\ 2 < \\v k — v l \\ 2 . Thus, u k converges in L 2 
to u*. 

The estimates we have already shown yield 

II^IIp < II^IIp < + 

for sufficiently large k. It means, that we can select a weakly convergent subsequence with limit 
u. Due to uniqueness of the limit, u* = u. Moreover, all weakly convergent subsequences have 
a common limit u*. Hence the sequence u k converges weakly in W 2 to u*. 

We know that J v has a unique minimizer u. Now, we have to show that u* is the minimizer 
of J v , i.e. u* = u. Obviously, we have 

J vk (u) > J v k(u k ). (5.9) 
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Due to the lower semicontinuity of the BV norm, we have 

rb \ f^ 1 

liminf J v k(u k ) = liminf h / \u k \+\im - / (u k -v k ) 2 > h / / (u*-v) 2 = J v (u*). 

k^oo k^oc J a fc-^oo 2 J a J a 2 J a 

On the other hand, we have 

lim J v k(u) = I h\u x \ + lim - / (u — v k ) 2 = J v (u). 

k-^oo J a fc^oo 2 J a 

Thus, due to (15.91) , we infer that 

J v (u) > J v {u*). 

Since u is a unique minimizer of J„, we conclude that u — u*. Our claim follows. □ 
We are ready to show our main results. 

Proof of Theorem \5.1\ Step 1. We have already noticed in Lemma I5TT1 that there exists a 
minimizer of Jh, v - Hence, there exists a solution to the following inclusion 

hdj{u) +u-v 3 0. 

Moreover, it is also unique, because if we had two, say u 1 and u 2 , then for some G dj(u l ), 
i = 1, 2, we had 

u 2 - u 1 + h(c 2 - C 1 ) = o. 

Once we apply the test function u 2 — u 1 to both sides, we see that \\u 2 — u 1 ^ < 0. Hence the 
claim, i.e. for any v E L 2 there exists u E D(J) a unique minimizer of J v . The above argument 
yields only that u belongs to BV. Now, the goal is to improve regularity of minimizers. 

Step 2. We will call by v e the standard mollification of v. Of course, ||t> e ||i jP < \\v \\i <p , but v e 
may not satisfy the boundary conditions, so we add a linear function. We call the result by v e . 
Of course, ||i> e ||i J p < ||^e||i,p + 0(e). 

We will show that the sequence of solutions u t to the minimization problem converges weakly 
in Wp and strongly in L 2 to u a solution to the original problem. 

Step 3. Since v e is smooth, then the sets E + (v e ), E~ (v e ) which we defined in Lemma [531 are 
open, i.e. 

E+(v e ) = |J (of, ft), E~(v e ) = |J (a~, (37). 

The index sets I + , I are at most countable. We may arrange the intervals at will. 

Step 4. We know by Lemma [531 above, that if v is smooth and the sets I + , I are finite, then 
u E Wp, for any p E [1, +oo) and it is piece-wise smooth. Moreover, v = uon E + (u) U E~(u). 
In particular, the set [a, b] \ E + (u) U E~(u) is a finite sum of closed intervals, so that we may 
write 

N 

[a,b}\E + (u)UE-(u) = \J[^,^]. 

i=i 

In particular, it is possible that £~ = £f. 
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We also know that if for some 5 > function v is monotone on [£r — 5, + 5], then u = v 
on £+], i.e. f ([^ r , * s a zero curvature facet. More interesting is the case, when for some 
5 > function v is convex or concave on [£f — 5, + 5]. Then, u = v (£,") = f (^ + ) on £+] 
and 

— v(x)) dx = 2h. 



From these properties, we deduce that 

\\ux\\p < \\v x \\p, (5.10) 

for all p G [1, oo). 

Step 5. In Lemma 1531 we constructed a sequence of continuous, piecewise smooth v k con- 
verging weakly to v e in W\. 

Let us call by u k the minimizers of J v k. Monotonicity of dj implies convergence of u k in 
L 2 . Indeed, if u k + dj(u k ) 3 then taking difference and applying it to a test vector yields, 

(Ue-U l e ,p) + h(j k -j l ,p) = (V k e -V l e ,p), 

where j k G dj(u k ), j l G dj(u l t ). When we choose p = u k — u{, then monotonicity of the 
subdifferential implies 

\\u k ~U% < ||v*-v'||2. 

Hence, the L 2 convergence of v k implies the L 2 convergence of u k to a limit u e . We have to 
improve the regularity of the limit. For this purpose, we notice that the estimate (15.101) applied 
to the sequence v k yields, 

\\ii k II < ll?; fc II 

|| a e,x\\P — II u e,x\\P 

for any p G (1,2]. Hence, we can select a weakly convergent subsequence in W\ with limit uf". 
Due to uniqueness of the limit we conclude that u e = uf", i.e. u e is in W\ for any finite p. This 
also implies that u k converges to u e uniformly. 

Since the norm is lower semicontinuous we also infer that 

1 1 fJ'e, x ||p — ||^e,a:||p — ||^x||p- 

So the same argument permits us to pass to the limit with e — > to conclude that u e converges to 
a limit u strongly in L 2 , and weakly in W 2 l - 

Step 6. We have to show that u e , for e > 0, and u are minimizers of J Vi for the corresponding 
data v € or v. For this purpose, we invoke Lemma I5T61 □ 

We also note a conclusion from the proof of Theorem 15. II 

Corollary 5.1. Let us suppose that v is continuous and piecewise smooth, such that one sided 
derivatives exit everywhere. The sets {v x > 0}, {v x < 0} are open with a finite number of 
connected components denoted by K. Then, u the unique minimizer of J v , belongs to Wp, for 
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any p G [1, +00) and it is piecewise smooth. Moreover, v = u on E + {u) U E (u) and there 
exists a G W}, such thata(x) G sgnu x (x) and 

h d 

—11— a = v — u. 

ax 

Furthermore, \\u x \\ p < \\v x \\ p , forallp G [1, 00). 

Theorem 15. H is slightly too general for our purposes, Theorem 15 .21 is its refinement. We will 
prove it momentarily. 

Proof of Theorem 15.21 Part (a) is obvious when K{v) = 00. If K(v) < 00, then the claim 
follows from the construction of u if h is sufficiently small. For a general h we have to use 
Lemma [54] 

Our proof of part (b) starts with the observation that v x G BV implies v x G L^. Hence, we 
can pass to the limit with p in the estimate || u \\ i )P < \\v\\x >p . Thus, \\u\\i + \\u x \\i < \\v\\x + \\v x \\i. 

If v x G BV, then by the general theory, see e.g. [0, there exists a sequence of smooth func- 
tions, {v k }, such that Ht^llsy converges to H^Hsv- We apply Lemma 1331 to deduce existence 
of a sequence {v km } such that the sets {v km , x > 0} and {v krri:X < 0} are open and have a finite 
number of components. Moreover, lim v km = v k in W\. 

Now, it is easy to calculate the norm Uiis-m^Hey f° r me corresponding minimizers u km for 
sufficiently small h. We have 

/ \Du km>x \ = [ \Dv km>x \ + J2(\ v km,Mt{h))\ + \v; miX {& +1 {h))\) 

< y~" / \Dv kniiX \ 

+ Y,^M( h y)- v L,MW)\ + Km,M+i(h)) - w(^+i(o))l) 

i 
i 

< E / i^^i + E(i v *™*& + (°))i + iv(^i(°))i) 

■ ■/« i + (0),& 1 (0)) 



/ \DVkm )X \. 
J a 



Here, we use the convention that if ^f(h) > a, then we write £o(h) = a and £^+1 = & provided 
that £ N < b. 

That is, we have 

|pfcro,a;||.BV < || u fcm,a: \\BV- 

We can find m k converging to zero as k goes to infinity such that Hf/cm^H w < ||^,a;||w + l/^- 
Finally, we use []Zl Theorem 5.2.1] to conclude that 

\\Du x \\ < liminf ||-Dw fcmfe;a; || < lim (HLh^H +l/k)= \\v x \\. □ 

k— ¥00 k— >oo 
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6 Asymptotics and examples 



Here, we present the proof of Theorem I2.3L an example of an explicit solution and numerical 
results describing the time behavior of solutions. 



6.1 A proof of Theorem 123 



Here is the argument. There is a finite number N of facet merging events 

= to < t\ < . . . < t/v < oo, 

when u has no time derivative but only the right-time derivative. Moreover, iV < K ess (uo, x )- We 
shall estimate maXj =0i ...jv-i{*i+i — U}. Let us set 

B = maxjaf,, a e }, b = minjaf,, a e }, Am = max-uo(^) — B, A m = b — mmuo(x), 

and £ = 1 is the length of / = [0, 1]. We notice that since our solution is almost classical, then u t 
exists except t E {t , t\, . . . , t^}. Moreover, u t is the vertical velocity of u. It is obvious from 
the definition of the composition 6 that the absolute value of (signow^.);,, is bigger or equal 2/1. 
We notice that the distance each essential facet travels in the vertical motion between collisions 
is no bigger than 

A = max{A M , A m , B - b}. 
Since we have a lower bound on the vertical velocity of u, we conclude that 

A 2 

max {t i+ i -U} < A - - 

i=0,..JV-l £ 

Thus, we have the following estimate 

text < 2K ess {u 0tX )A/£. (6.1) 

Hence K ess (u x (t ext ) = 0, then thus u{t) for t > t ext is a monotone function being a stationary 
state of the system. □ 

6.2 An explicit solution 

In order to illustrate the behavior of a particular solution we take x 2 as an initial datum for (11.11) . 
We consider this system on the interval ( — 1,1), 

u t -£sgnu x = in (-1, 1) x (0, T), 

u(-l,t) = u(l,t) = 1 for te(0,T), (6.2) 

u\ t=0 = x 2 for (-M)- 

The proved results quarantee us the following form of the solution to (16.21) . 

, , / a 2 (t) for \x\ < a(t), ( 

UVX,t) — { 2 r 11/ f±\ i (6-3) 

v ' \ x 2 for \x\ e (a(t), 1) v 
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By Definition [24] we get that 



^sgno^| ha(tMt)] = -L 



Thus by (16.21) and (16.31) we find a relation on a(t) as follows 

1 l~3~ 
d t a 2 (t) = hence a(t) = y -t 

to keep the agreement to the initial datum. 

Summing up the length of the facet is 2a{t) = 2y|t, the speed of it is d t a(t) ~ t~ 2 l z and 
the extinction time of u = 1 is T sta b = f • 



6.3 Numerical simulations 

Now, we are prepared to computer implementations of our results. Simulations were done in 
Octave package. The main part of the program is a loop running until the graph reaches it's final 
shape. During one step all facets (i.e. points where E df) are moved until (if it is possible) 
each of them fills the area equal to 2h. In the pictures shown below we used h = 5. The reason 
why it may be not possible to fill the 2h area is that the moving facet may reach the boundary of 
the interval that it is defined on or it may reach the boundary of another facet after it filled the 
required area (whereas each of them moved separately may fit its domain). When any of these 
interactions happens, we change the h value for a maximum reached value (let us call this new 
value hmin) and move all facets so that they fill the area of 2h min . We use h min just in this one 
step but for all facets and then get back to h value. After each step, we recalculate domains and 
check if we still use all functions (some of them may disappear, as the x 2 — 2x function defined 
on [0, 1] interval after the first step of the v\ example from table 1). 

In none of the presented examples a facet fills the maximum area. We chose h big enough to 
avoid unnecessary steps. 

We calculate the time a step takes as 2h ^ n ■ We do this using the following logic — we make 
an assumption that one full step (i.e. area of 2h is filled) is my time unit, two full steps count as 
t = 2, ~/t takes t = | to fill. In the pictures accumulated time is presented. 

As an initial data in three presented examples, we use functions described in the table below. 
The first column contains intervals which set the domain, the next three columns contain formulas 
for respective examples: 
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domain 




v 2 


f ^ 


[-1.5,-1] 


x-2 


3x 2 + llx 


3x 2 + llx 


[-1,0] 


-x 2 + x + 2 


—x 2 + 5x + 1 


—x 2 + 5x + 1 


[0,1] 


x 2 — 2x 


x-2 





[1,2] 


-x 2 + 5x + 1 


2x-7 


2x- 7 


[2,3] 


x 2 — Qx + 8 


x 2 — Qx + 8 


1 


[3,4] 





-x 2 + x + 2 


-x 2 + x + 2 


[4,5] 


2x - 7 


x 2 — 2x 


x-2 


[5,5.5] 


1 


x 2 + 15x 


x 2 + 15x 



Table 1. Examples 1, 2, 3 (respectively) used in the simulations 



To create the three examples, we use the same domain and permute functions to obtain in- 
teresting shape. In some cases, we have to move parts defined on some intervals vertically to 
obtain continuous result. Therefore, in some cases the same function used on the same interval 
has different values. What is more, we move the whole graph vertically so that the smallest value 
is 1; it makes integration easier without changing the shape of solutions. 

We use polynomials as an approximation of a continuous function defined on closed inter- 
val; in the examples mentioned they are of degree 2, but the algorithm remains the same for 
polynomials of higher degree. Functions defined on intervals model situation of non-continuous 
derivative. 

Let us look at results of simulations presented on the figures: 




-1 1 2 3 4 5 6 



Figure 1 : The first example 
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t=0.00 



t=0.06 



t=0.24 



t=0.33 




-1 1 2 3 4 5 



Figure 2: The second example 




-1 1 2 3 4 5 



Figure 3: The third example 
Observe that, all degenerated facets disappear after the first step of evolution. The number of 
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regular facets that may appear is limited by their number and the overall number of regular facets 
decreases from the second step of evolution. The flat area broadens with each step. All solutions 
remain continuous and their 1 1 ■ 1 1^ norm is bounded by the norm of initial data. 
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